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Abstract. In this paper we give, for the first time, a qualitative description of the 
asymptotic dynamics of a class of non-tilted spatially homogeneous (SH) cosmologies, 
the so-called exceptional Bianchi cosmologies, which are of Bianchi type VI_i/ 9 . This 
class is of interest for two reasons. Firstly, it is generic within the class of non-tilted SH 
cosmologies, being of the same generality as the models of Bianchi types VIII and IX. 
Secondly, it is the SH limit of a generic class of spatially inhomogeneous G 2 cosmologies. 

Using the orthonormal frame formalism and Hubble-normalized variables, we show 
that the exceptional Bianchi cosmologies differ from the non-exceptional Bianchi 
cosmologies of type Vlh in two significant ways. Firstly, the models exhibit an 
oscillatory approach to the initial singularity and hence are not asymptotically self¬ 
similar. Secondly, at late times, although the models are asymptotically self-similar, the 
future attractor for the vacuum-dominated models is the so-called Robinson-Trautman 
SH model instead of the vacuum SH plane wave models. 


PACS numbers: 0420H, 0440N, 9880H 

1. Introduction 

In this paper we give a qualitative description of the asymptotic dynamics of a class 
of spatially homogeneous (SH) cosmological models which to date have received little 
attention in the literature. The matter content is a perfect fluid with equation of state 
p = (7—l)/r, where 7 is a constant. The models are non-tilted, i.e. the fluid 4-velocity is 
orthogonal to the orbits of the 3-parameter group G3 of isometries, that is of Bianchi 
type VI/j, with the value of the group parameter h = This exceptional value 

permits the shear tensor to have an additional independent component (see Ellis and 
MacCallum 1969, p 123, case Bbii) with the result that this class of models depends 
on flve arbitrary parameters instead of four, making it of maximum generality, i.e. of 
the same generality as the non-tilted Bianchi VIII and Bianchi IX models. We shall 
refer to this class of SH cosmologies as the exceptional Bianchi cosmologies. 

The longer term goal in studying the dynamics of SH cosmologies is to gain insight 
into the dynamics of spatially inhomogeneous cosmologies. Indeed our interest in the 
exceptional Bianchi cosmologies stems from the fact that they are closely related to the 
simplest class of spatially inhomogeneous cosmologies, the so-called G2 cosmologies, 
which admit a 2-parameter Abelian isometry group G2. The relation is due to the fact 
that the present class^ admits an Abelian G2 of isometries as a subgroup of the full 

t Models of all Bianchi types except for types VHI and IX have this property. 
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isometry group. Most research on G 2 cosmologies has dealt with two special cases, 
either diagonal G '2 cosmologies, in which case the line element is diagonal, or more 
generally, orthogonally transitive G 2 cosmologies, in which case the line element has a 
2x2 block form (see van Elst et al 2002, where other references can be found). The 
exceptional Bianchi cosmologies are unique among the non-tilted SH cosmologies in 
that the G 2 subgroup does not act orthogonally transitively, and, as a result, this class 
of models is the SH limit|| of the most general G 2 cosmologies. One thus expects that 
the asymptotic dynamics of the present class of SH cosmologies will provide insight 
into the dynamics of general G 2 cosmologies. 

In this paper we show that the presence of the additional shear degree of 
freedom has a significant effect on the asymptotic dynamics of the exceptional Bianchi 
cosmologies, both in the singular regime and in the late-time regime. Firstly, this shear 
degree of freedom destabilizes the vacuum Kasner solutions, leading to an oscillatory 
singularity. Secondly, it destabilizes the SH vacuum plane wave solutions, with the 
result that they no longer play a dominant role as regards the late-time dynamics. 

The plan of this paper is as follows. In section || we present the evolution equations 
for the exceptional Bianchi cosmologies, using the orthonormal frame formalism, and 
discuss some properties of the Hubble-normalized state space. Details of the derivation 
are provided in appendix A. In section ^we give the equilibrium points of the evolution 
equations, and describe their local stability, indicating the various bifurcations that 
are determined by the equation of state parameter 7 . Sections ^ and ^ contain the 
main results of the paper, a description of the asymptotic dynamics of the models at 
late times and near the initial singularity, in terms of the future attractor and the past 
attractor of the evolution equations. Finally, in section || we conclude by pointing out 
some analogies between the exceptional Bianchi cosmologies and other classes of SH 
cosmologies with an oscillatory singularity. We also comment on the implications of 
our results for G '2 cosmologies. 

It is assumed that the reader is familiar both with the orthonormal frame 
formalism of Ellis and MacCallum (see MacCallum 1973), and with the use of 
dynamical systems techniques in cosmology, for which we refer the reader to the book 
of Wainwright and Ellis (1997)|. We use geometrized units with SttG = c = 1 and the 
sign conventions of MacCallum (1973). 

2. Evolution equations 

In this section we give the evolution equations for the exceptional Bianchi cosmologies. 
The equations are derived in appendix A using the orthonormal frame formalism and 
Hubble-normalized variables. The models are described by a dimensionless state vector 

a; = (E+,E_,S2,Sx,IV_,A), (2.1) 

subject to a constraint given by 

g{x) = (E++ VST,-)A-T,^N-=0. ( 2 . 2 ) 

The variables E-|_, S_, E 2 and Ex describe the shear of the fluid congruence, while 
N- and A describe the spatial curvature. These variables are dimensionless, having 

f There are also classes of tilted SH models that are SH limits of the most general G 2 cosmologies 
(see Wainwright and Ellis 1997, table 12.4, p 268). 

I Henceforth, we will refer to this work as WE. 
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been normalized with the Hubble scalar^ H, which is related to the overall length 
scale £ by 


H = 


(2.3) 


The overdot denotes differentiation with respect to clock time along the fluid 
congruence. The state variables depend on a dimensionless time variable r that is 
related to the length scale £ by 

£ = £oe^, (2.4) 

where £o is a constant. The dimensionless time r is related to the clock time t by 
dt 1 
d)^ ^ F’ 


(2.5) 


as follows from equations (2.3) and (2.4). In formulating the evolution equations we 
require the deceleration parameter q, defined by 


and the density parameter H, defined by 


n = 


3i?2- 


( 2 . 6 ) 


(2.7) 


The evolution equations for the components of a; in (^) take the following form 

E+ = (g - 2)S+ + 3E^ - 2Nf - 

E'_ = (g- 2)S_ - V5El+2VST,l - 2V3Nf + 2V3 , 

= (g-3E+ + V3E_-2)E2, 

E; = (g- 2V3E_ - 2)Ex - 8N-A, 

N'_ = (g + 2E+ + 2v^E_)Af_ +6ExH 
A' = (g + 2E+)H, 


( 2 . 8 ) 


where 


g =2E2 + i(37-2)H, 

Vl =l-Tf^ -Nf- 4A^, 


E^ = E^ 




(2.9) 
( 2 . 10 ) 
( 2 . 11 ) 

and ' denotes differentiation with respect to r. There are two auxiliary equations, the 
evolution equation for H, which reads 

= [2g - (37 - 2)]H, (2.12) 

and the evolution equation for the constraint function g(a;), namely, 

g' = 2(g + E+-l)g, (2.13) 

which guarantees that g(x) = 0 defines an invariant set. Both (2. il) and ( ^.13| ) follow 
from (2.8). 

t On account of (2.3), if is related to the rate of volume expansion © of the fluid congruence according 

to H= ie. 
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The physical state space is the subset of ' 
the requirement 


defined by the constraint (2.2) and 


^ 0. (2.14) 

The restriction ( p.l4 ) implies that the state space is bounded. In addition, the 
evolution equations are invariant under the symmetries 


(S+,S_,S2,Sx,iV-,4l) ^ (S+,E_,±S2,±Sx,±iV_,±A), 


(2.15) 


provided that the product T,xN-A does not change sign. Since the evolution equations 
for Yj 2 and A do not allow either of these variables to change sign along an orbit, we 
can assume without loss of generality that 


E2 ^ 0, A ^ 0. 


(2.16) 


Thus the physical state space T) is the subset I? C K.® defined 

(Hi)- 


by (lU), (^-141 ) and 


The state space T) contains a number of invariant subsets which play an important 
role as regards the dynamics in the asymptotic regimes. Firstly there is the vacuum 
boundary, given by = 0. Secondly there are various invariant sets with A = 0 that 
describe vacuum Bianchi I (Kasner) and vacuum Bianchi II (Taub) models, which we 
shall introduce in section ||. Thirdly, it follows from the evolution equations ( p.Sj ) that 
the condition^ 


Ex=0, N_=0, E+ + v^E_ = 0, 


(2.17) 


define a 3-dimensional invariant set, which we shall denote by S. This set describes 
the subclass of the exceptional Bianchi cosmologies for which the G 2 admits one 
hypersurface-orthogonal Killing vector held. The asymptotic dynamics of this subclass 
have been analyzed in detail by Hewitt (1991). Finally, we note that the invariant 
set E 2 = 0 describes the non-exceptional Bianchi VI_i /9 models. In other words, 
the variable E 2 describes the additional shear degree of freedom referred to in the 
introduction. 


3. Local stability of the equilibrium points 


In this section we discuss the local stability of the equilibrium points of the DE (2.8) 
subject to the constraint (^.2|). The equilibrium points can be found by solving the 
system of algebraic equations 


f{a) = 0 , g{a) = 0. 

The local stability is determined by linearizing the DE (|^) x — a, which gives 
x' = T)f{a)x, 


and hnding the eigenvalues of the derivative matrix D/(a). The existence of the 
constraint (2.2) complicates the analysis, which requires that we consider only the 
physical eigenvectors of D/(a), that is, those which are tangent to the constraint 
surface, or equivalently, those which are orthogonal to the gradient vector Vg{a). 
Eigenvalues and eigenvectors that satisfy this condition shall be referred to as physical. 
We note that if all the physical eigenvalues have negative (positive) real parts then 
the equilibrium point is a local sink (source), that is, it attracts (repels) all orbits in 


f On account of the constraint ( |2.2|) the first two restrictions imply that + \/3E_ = 0, since 
A A 0 for the exceptional Bianchi cosmologies. 














Asymptotic dynamics of the exceptional Bianchi cosmologies 


5 


a neighbourhood. In addition to isolated equilibrium points, we will also encounter 
arcs of equilibrium points, for which one physical eigenvalue is necessarily zero. In 
this case, the criterion for a local sink (source) is that all eigenvalues other than the 
zero one have negative (positive) real parts. 

We now list the equilibrium points of the DE (2.8) subject to the constraint (2.2), 
and give the values of the density parameter and the deceleration parameter q. We 
refer to WE (see pp 189-93) for the corresponding exact solutions. 


Non-vacuum equilibrium points (Tl > OJ 

3.1. Flat Friedmann-Lemaitre equilibrium point, FL 


E-)- = E_ = E 2 = Ex = N— = A = 0, 
n = l, g=i(37-2), 0<7<2. 


3.2. Collins-Stewart Bianchi II equilibrium points, CS^ 

E+ = -^(37-2), E_ = V3E+, iV_ = ±1^73(37-2)(2- 7 ), 

E 2 = Ex = A = 0, 

^ = 9=5(37-2), I <7 <2. 

3.3. Collins Bianchi V 7 _i /9 equilibrium point, C 


E+ = -i(37-2), S_=-^E+, v4= 1^7(37-2)(2- 7 ), 

E 2 = Ex = IV- = 0, 

fI = i(5-37), g=i(37-2), | < 7 < |. 


3 . 4 . Wainwright Bianchi FJ-i/g arc of equilibrium points, W 


S+ 


1 

3’ 



S 2 


\/5 a 

3 V 3 ’ 


Ex = = 0, A= ^i(4 + 5a2), 0 < a < 1, 

fl = |(l-a2), q = l, j=f. 


Vacuum equilibrium points (fl = 0) 

3.5. Robinson-Trautman Bianchi V7_i/g equilibrium point, RT 





S 2 = 


3V3’ 
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Table 1. Local sinks in the exceptional Bianchi state space. 


Range of 7 

Local sink 

Q 


w 

0<7=S § 

FL 

1 

0 

0 

3 ^ ^ ^ 9 

C 

4(5-37) 

1t(37-2)2 

4(37-2)v4TTT 


W 


^(4 + 5a^) 

^y/(4 + 5a-^)(13 + 35a'^) 

f <7<2 

RT 

0 

1 

3 

4 

3 


3.6. Plane wave arcs of equilibrium points, PW^ 

E+ = -i(4-a), S_ = iV3a, = 0, 

Ex = ±i-\/a(l - a), iV_ = -Ex, A = ^a, 0 < a < 1, 

q = 5 ( 4 - a). 

3. 7. Kasner circle of equilibrium points, JC 

E 2 = Ex = = A = 0 , E^ + E?. = 1 , 

9 = 2. 

We now consider whether any of the equilibrium points are local sinks or sources. 
It turns out that for each value of the equation of state parameter 7 in the interval 
0 < 7 < 2 , there is a unique equilibrium point/set that is a local sink, as indicated in 
table 1^. We can describe the transitions, i.e. bifurcations, that occur between these 
local sinks as follows: 

-7=2 7=10 7=10 

FL ^ C —? W -4 RT (3.1) 

The mechanisms for these bifurcations, without giving full details, can be 
described as follows. The linearization of the evolution equations for iV_ and A at the 
flat FL point is 

W = i(37-2)fV_, A'= 1(37-2)4, 

establishing that the spatial curvature variables N- and A destabilize FL at 7 = |. 
Moreover when 7 = |, both the Collins-Stewart CS and Collins C equilibrium points 
bifurcate from the FL equilibrium point. The linearization of the evolution equation 
for E 2 at the Collins equilibrium point C is 

E' =i(97-10)E2, 

showing that E 2 destabilizes C at 7 = Indeed, at 7 = ^ there is a line bifurcation 
which is characterized by the exchange of stability between the Collins equilibrium 
point C and the Robinson-Trautman vacuum equilibrium point RT by means of the 
Wainwright arc of equilibria W which connects both points. It is also of interest to 
note that for saddle, and it merges with the plane wave arcs (at 

a = I) when 7 = §■ 

The local sinks in the exceptional Bianchi state space are listed in table |^, together 
with the values of the density parameter fl, the shear parameter E, which describes the 
anisotropy in the fluid congruence (see WE, p 114) and the Weyl curvature parameter 
W, which can be regarded as a measure of the intrinsic anisotropy of the gravitational 
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field (see Wainwright et al 1999, p 2580). For the present class of models, and 
are given by ( 2.10| ) and (2.11). The formula for W is provided in appendix B. 

It is of interest to note that the plane wave arcs of equilibrium points, PW=*=, 
which describe the SH plane wave vacuum solutions, do not appear in the table. This 
result represents an important difference between the exceptional and non-exceptional 
Bianchi VI/i models. For the latter class, the plane wave equilibrium points are a local 
sink for vacuum models and for perfect fluid models, subject to a restriction on the 
equation of state parameter 7 . On the other hand, the plane wave arcs of equilibria are 
saddles in the exceptional Bianchi state space. To see this, we linearize the evolution 
equation for 112 at the plane wave arcs obtaining 


S 2 = ^6 — a)T,2- 

Thus, the extra shear degree of freedom, S 2 , that is present in the exceptional Bianchi 
cosmologies, destabilizes the plane wave solutions. Their stability is, in effect, inherited 
by the Robinson-Trautman solution. 

As regards local sources, it follows upon analysis of the eigenvalues associated 
with the equilibrium points that there are no local sources in the exceptional Bianchi 
state space. In particular, the Kasner circle /C is a saddle, and, consequently, the initial 
state of a typical model cannot be a single Kasner equilibrium point. 


4. The late-time asymptotic regime 


We have seen that for each value of 7 in the range 0 < 7 < 2, excluding 7 = there 
is a unique equilibrium point that is a local sink of the evolution equations, while if 
7 = ^, the local sink is an arc of equilibrium points. We note that all of the local 
sinks in table lie in the three-dimensional invariant set S defined by (2.17). It has 
been proved that each of these local sinks is the future attractor in the invariant set 
S, for the relevant range of 7 (see Hewitt 1991)|. 

We have performed extensive numerical experiments which suggest that the 
solution determined by an arbitrary initial condition satisfies 

lim Sx = 0, lim N- = 0, lim E_|_ -|- -s/SE- = 0. (4.1) 

T^ + CX) T^ + OO T—* + 00 


In other words, the experiments provide strong evidence that all orbits are attracted 
to the invariant set S, i.e. that the future attractor is contained in this subset. We 
thus make the following conjecture. 

Conjecture 4.1. Each local sink given in table |^/or 7 in the range 0 < 7 < 2 is the 
future attractor in the exceptional Bianchi state space. 

This conjecture is in fact known to be true for 7 in the range 0 < 7 < | (see WE, 
theorem 8.2, p 174). Giv ing a proof for 7 in the range | ^ 7 < 2 necessitates 
establishing the limits ( 0 ). 


5. The singular asymptotic regime 

In this section we show that the exceptional Bianchi cosmologies exhibit an oscillatory 
approach to the initial singularity, as do the Bianchi VIII and IX models. There are 
two significant differences, however, which we now describe. 

f The analysis in this paper is complete, except for the case 7 = ^ ■ We have recently been able to 
complete this analysis by finding a Dulac function for the family of two-dimensional invariant subsets 
that foliate the state space. We thank Neville Dubash for assistance with this matter. 
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Figure 1. The arrays + H-, etc. give the signs of the eigenvalues Xff_ , ^E 2 ^.nd 

in that order. The variables(s) listed next to each of the six arcs indicates which of the 
variable(s) N—, S 2 , Sx is growing into the past. Kasner points on the bold arcs T2Q1 
and T3Q2 have a two-dimensional Kasner transition set into the past. 


In the dynamical systems approach, the mechanism for creating an oscillatory 
singularity is that the Kasner equilibrium points, which comprise the Kasner circle, are 
saddles, and the unstable manifold (into the past) of any Kasner point is asymptotic 
to another Kasner point. In other words, the unstable manifold consists of orbits that 
join two Kasner points. We shall refer to these unstable manifolds as Kasner transition 
sets, because they provide a mechanism for a cosmological model to make a transition 
from one (approximate) Kasner state to another, as it evolves into the past. For non- 
tilted Bianchi VIII and IX models, each Kasner transition set is one-dimensional and 
thus there is a unique heteroclinic orbit|| joining two Kasner points. In the present 
case, for two of the six arcs on the Kasner circle, the Kasner transition set is two- 
dimensional, leading to a one-parameter family of heteroclinic orbits. Nevertheless, 
in both cases, the stability properties of the Kasner equilibrium points lead to the 
existence of infinite heteroclinic sequences, i.e. infinite sequences of equilibrium points 
on the Kasner circle joined by special heteroclinic orbits, oriented into the past (see 
WE, section 6.4.2, for Bianchi VIII and IX models). These heteroclinic sequences 
govern the dynamics in the singular regime (t — 00 ) in the sense that a typical 
orbit shadows (i.e. is approximated by) a heteroclinic sequence as r ^ — 00 . In 
physical terms, the corresponding cosmological model is approximated by a sequence 
of Kasner vacuum models as the singularity is approached into the past, the so-called 
Mixmaster oscillatory singularity. 

The second difference is that in the case of Bianchi VIII and IX models, the 
instability of the Kasner points is generated by spatial curvature only, while in the 
case of exceptional Bianchi cosmologies, the instability is due to a combination of 
spatial curvature and off-diagonal shear, in particular, the variables N-, E 2 and Ex. 

Figure shows the signs of the eigenvalues of the Kasner equilibrium points 
associated with the variables V_, E 2 and Ex, showing that points on the arcs 
T 2 Q 1 and T 3 Q 2 have a two-dimensional Kasner transition set into the past (negative 
eigenvalues indicate instability into the past). The figure also shows which of the three 
variables are increasing into the past in a neighbourhood of the Kasner circle. 

t A heteroclinic orbit is one that joins two equilibrium points. 
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z- l^:. iz. 



past. 


Figure|| shows the projections in the E+I]_-plane, of three families of heteroclinic 
orbits, that join two Kasner points. For a Kasner point with a one-dimensional 
transition set, the orbit shown in figure ^ emanating from that point is the transition 
set. These families are described as follows. 

Sn_-' Vacuum Bianchi II models (Taub models) 

These orbits describe the familiar vacuum Bianchi II models. They are given by 

E 2 = Ex=A = 0, ^ iV-7^0, (5.1) 

with E_ -I- -s/S = to(E+ -|- 1), where to > -^ is a constant. 

Ssyf • Kasner models relative to a non-Fermi-propagated frame, with E 2 = 0 
These orbits are given by 

^^=N_=A = 0, III + = 1, ^xi^O, ( 5 . 2 ) 

with E+ = TO, where —l<TO<lisa constant. 

.' Kasner models relative to a non-Fermi-propagated frame, with E x = 0 
These orbits are given by 

Ex=iV_=A = 0, El + El+^l = l, E2>0, (5.3) 

with E+ -|- E_ = TO, where —2 < to < 2 is a constant. 

Figure^ gives a representation of the two-dimensional transition set of a Kasner 
point on the arc T^Q 2 - This invariant set consists of a one-parameter family of orbits 
with n = 0, A = 0, Ex =0, N- 7^ 0 and Ex 7^ 0, which join Aio B (the projection of 
two typical orbits is shown), and two special orbits, one with N- 7 ^ 0 and E 2 = Ex =0 
(a member of the family Sn_), and the other with E 2 7 ^ 0 and N- = Ex = 0 
(a member of the family S^^)- These orbits, which describe the vacuum Bianchi II 
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E. 



Figure 3. The two-dimensional Kasner transi¬ 
tion set of a Kasner point A on the arc T 3 Q 2 , pro¬ 
jected onto the E+S_-plane. The arrows show 
evolution into the past. 


E. 



Figure 4. The two-dimensional Kasner transi¬ 
tion set of a Kasner point A on the arc T 2 Q 1 , pro¬ 
jected onto the S+E_-plane. The arrows show 
evolution into the past. 


models relative to a non-Fermi-propagated frame, can be described explicitly by the 
equation|f| 

(1 + S+ + V3E_)2 + 3V2 

- 7 = -= 0^1 (5.4) 

(4 + E+ + V3E_)2 

where a is a constant that depends on the point A. 

Figure^ gives a representation of the two-dimensional transition set of a Kasner 
point on the arc T 2 Q 1 . This invariant set consists of a one-parameter family of orbits 
with n = 0, A = 0, N- = 0, E 2 ^ 0 and Sx 0, which join Aio B (the projection of 
two typical orbits is shown), and two special orbits, one with E 2 ^ 0 and Sx = iV_ = 0 
(a member of the family ), and the other with S x 7 ^ 0 and N_ = S 2 = 0 (a member 
of the family 5s All orbits in this manifold, which intersect the Kasner circle in 
six points, as shown in hgure represent one and the same physical Kasner model. 
The six Kasner points referred to have the same Kasner exponents, but differ in the 
labelling, and the orbits in the Kasner transition set thus describe a rotation of the 
spatial axes relative to which this Kasner model is described. These orbits can be 
described explicitly by the equation|| 

2 - 6 E+ -h 8S3 -p 373(^3 E+ -h S.jSf = a, (5.5) 

where a is a constant that depends on the point A, thereby distinguishing different 
Kasner models. 

Numerical experiments suggest that as r ^ — oo, an orbit that is close to a point 
A on the arc T 3 Q 2 will shadow either the orbit AD C Sn_ or the orbit AC C (Sss, 

f This expression is a first integral of the evolution equations ( [2.^ ) when restricted to the invariant 
set Q = 0, A = 0, Sx =0. An equivalent expression has arisen in the analysis of G 2 cosmologies (see 
Weaver 2002). We thank Marsha Weaver for poi ntin g out this first integral in the G 2 context, and 
Woei Chet Lim for calculating the expression in ( p.4j ). 

t This expression is a first integral of the evolution equations ( [2.^ ) when restricted to the invariant 
set Q = 0, A = 0, N— = 0. It is related to a Hubble-normalized scalar formed from the Weyl tensor, 
which is constant for Kasner models. We thank Woei Chet Lim for verifying this result. 
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rather than shadowing a typical orbit in the two-dimensional Kasner transition set. 
A similar remark applies to points on the arc T 2 Qi. The result of this behaviour is 
that as T —oo an orbit will shadow a heteroclinic sequence that consists of orbits 
belonging to the invariant sets Sn_, and This behaviour motivates the 

following conjecture concerning the past attractor A~ . 

Conjecture 5.1. The past attractor is the two-dimensional invariant set consisting 
of all orbits in the invariant sets Sm_, and figure^ and the Kasner 

equilibrium points, i.e. 

A = S]\f_ U Sy:,, U Sy2 U /C. (5.6) 


This conjecture can be formulated in terms of limits of the state variables as 
follows. Referring to (|2.10D and ( |2.1lD we see that the set A is defined 

by 

rj = 0, A = o 


and 


N_Y.2 = 0 , 


S 2 SX — 0, 


SxAr_ =0, 


the final equality being a consequence of A = 0 and the constraint (2.2). Thus our 
conjecture concerning the past attractor can be formulated as 

lim 12 = 0 , 


lim A = 0, 

r—»• —00 


where 


A = (A_I] 2 )" + (S 2 Sxf + (SxA^-)". 


These limits imply that limT^-oo A = 0, on account of (^). Note that for a generic 
orbit, lim 7 .^_oo(A^-, S 2 , Sx) does not exist. 


6. Discussion 

There are two generic classes of non-tilted ever-expanding SH cosmologies, namely 

i) the Bianchi type VIII cosmologies, and 

ii) the exceptional Bianchi cosmologies. 

The evolution of the first class in the asymptotic regimes has been described in 
the literature (see for example, WE, pp 143-7 for the singular regime, and the 
accompanying paper, Horwood et al 2002, for the late-time regime). In this paper 
we have used Hubble-normalized variables to describe the asymptotic dynamics of the 
second class. We have shown that with this choice of variables, the Einstein field 
equations reduce to an autonomous DE on a 5-dimensional compact subset of K.®. 
Since the state space is compact, there is a past attractor and a future attractor, which 
we have been able to identify using heuristic arguments and numerical experiments, 
thereby describing the dynamics in the singular and late-time asymptotic regimes. 

The mathematical description of the exceptional Bianchi cosmologies differs from 
that of the non-exceptional models in the following ways: 

i) the expansion-normalized state space is of dimension five rather than four, 
corresponding to the presence of an additional shear degree of freedom, 

ii) the Abelian G 2 subgroup of the isometry group does not act orthogonally 
transitively. 
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We have shown that these differences lead to significant changes in the asymptotic 
dynamics: 

i) the singular regime is oscillatory (infinite sequence of Kasner states), rather than 
asymptotically self-similar, 

ii) there is a new asymptotic state described by the vacuum SH Robinson-Trautman 
solution, which describes the late-time behaviour of vacuum models and of perfect 
fluid models with 7 satisfying ^ < 7 < 2 . 

These two points also highlight the similarities and the differences between the two 
generic classes of ever-expanding models: both classes have an oscillatory singularity, 
but differ as regards the late-time asymptotic regime. The Bianchi VIII models are 
not asymptotically self-similar at late times and exhibit Weyl curvature dominance 
(the Weyl curvature parameter W diverges as t ^ -boo, see Horwood et al 2002), 
whereas the exceptional Bianchi cosmologies are asymptotically self-similar and hence 
do not exhibit Weyl curvature dominance (W has a finite limit as r ^ -boo, given in 
table 1 ^. 

We now discuss some aspects of the oscillatory singular regime that occurs in 
SH cosmologies. This type of singularity, also known as a Mixmaster singularity, was 
first discovered in vacuum Bianchi IX cosmologies (see Misner 1969 and Belinskii et 
al 1970) and subsequently in the closely related Bianchi VIII cosmologies (see for 
example, WE, section 6.4). Within the dynamical systems framework, the oscillatory 
singularity in these models is created by the three degrees of freedom in the spatial 
curvature, which destabilize the Kasner equilibrium points in the Hubble-normalized 
state space. In the present models, there is only one degree of freedom in the spatial 
curvature, but instead there appear two off-diagonal shear degrees of freedom, which 
again destabilize the Kasner equilibrium points. This destabilization, and the resulting 
existence of the Kasner transition sets (see section ||) , leads to the creation of infinite 
sequences of heteroclinic orbits joining Kasner equilibrium points, resulting in an 
oscillatory singularity. This phenomenon can in fact occur in models of all Bianchi 
types, if one permits sufficient general source terms for the gravitational field. Two 
classes which have been analyzed in detail and are closely related to the present class, 
are the Bianchi I magnetic cosmologies (see LeBlanc 1997) which have one magnetic 
degree of freedom and two off-diagonal shear degrees of freedom, and the tilted perfect 
fluid Bianchi H cosmologies (see Hewitt et al 2001) which have one spatial curvature 
degree of freedom and two off-diagonal shear degrees of freedom. Indeed, the common 
feature of all SH models with an oscillatory singularity is the occurrence of the above 
Kasner destabilization, which can be caused by any combination of spatial curvature, 
off-diagonal shear or magnetic degrees of freedom totaling at least three. 

We conclude by discussing the implications of our results for G 2 cosmologies. 
As mentioned in the introduction, the present class of SH cosmologies are the SH 
limit of the most general G 2 cosmologies, namely those for which the G 2 does not act 
orthogonally transitively. One thus expects that the initial singularity in a general 
G 2 model will be oscillatory, and indeed numerical evidence has been provided that 
suggests this is the case, at least for vacuum models (see Berger et al 2001 and 
Lim 2002). Secondly, the fact that the SH Robinson-Trautman solution is a local 
sink in the Hubble-normalized state space for vacuum models and for perfect fluid 
models with equation of state parameter 7 satisfying ^ < 7 < 2 suggests that this 
solution may play a role in the dynamics of general G 2 cosmologies at late times. It 
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would thus be of interest to investigate the stability of the SH Robinson-Trautman 
solution in this context. 
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Appendix A. Derivation of the evolution equations 


In this appendix we write the EFEs for the exceptional Bianchi cosmologies as 
an autonomous system of ordinary differential equations subject to one algebraic 
constraint. Our starting point is the system of evolution equations for a general 
spatially homogeneous cosmological model given in Hewitt et al (2001) (see equations 
(A.11)-(A.18)). These equations have been obtained by writing the orthonormal frame 
equations for a spatially homogeneous model using a group invariant orthonormal 
frame {eo,eQ,}, choosing Cq to be the unit normal to the group orbits. Furthermore, 
all commutation functions and source terms are normalized with the Hubble scalar H. 

We now describe the specializations that lead to the exceptional Bianchi 
cosmologies. Firstly, the models are non-tilted in the sense that the unit normal 
eo is aligned with the fluid velocity u. Furthermore, the models have a perfect fluid 
source with a 7 -law equation of state given by 

P = (7-1)U, 


where P is the Hubble-normalized pressure and H is the dimensionless density 
parameter. The equation of state parameter 7 is restricted to the interval 0 < 7 < 2. 
Secondly, the orthonormal frame vectors 62 and 63 are chosen to be tangential to the 
orbits of the two-dimensional Abelian subgroup G 2 of the three-dimensional symmetry 
group G 3 . The consequences of this frame choice are 

A^la = 0 , Aq, = (5^q,Ai, i ?3 = —S 12 , i ?2 = S 13 (A.l) 


(see Wainwright 1979, pp 2019-21). Equation (A.18) in Hewitt et al (2001) contains 
three constraints, one of which will give rise to the constraint ( 2 . 2 ) and the other two 
which together lead to the following possibilities: 


i) S12 — S13 — 0, 
or 

ii) 9 A 3 — A 23 -f N 22 N 33 = 0 and A^ 22 Si 2 + {N 23 — 3 Ai)Ei 3 = 0 . 

Making the former choice is equivalent to insisting that the Abelian G 2 subgroup acts 
orthogonally transitively. We make the latter choice, which results in the models being 
Bianchi VR with^ h = 

t Recall that the group parameter h is related to Aa a.nd according to 
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The remaining gauge freedom at this point is an initial alignment in the G 2 
orbits. We choose to set E 12 = 0 initially. It follows that this choice is preserved by 
the evolution equation for E 12 provided that 

i?i=E23. (A.2) 

The gauge freedom in the angular velocity variables Ra of the spatial frame allows us 
to make such a choice for i?i. Thus, the constraints in case ii) now yield 

N 23 - 3Ai = 0, N 33 = 0. (A.3) 

At this stage, the independent non-zero Hubble-normalized variables are 

S22, S33, S13, E23, N22, Ai. 


We relabel the shear variables according to 

S+ = 5(5122 + 5I33), S- = ^;^(S22 — 5I33; 

and use the notation 

N_ = ^ N 2 : 


5I2 — 51i3, Ex — —7^ E23, 


1 

V3 ■ 


'' 22 , 


A = Ai. 


(A.4) 


(A.5) 


The evolution equations (2.8) and the constraint (2.2) for the exceptional Bianchi 
cosmologies now follow from equations (A.10)-(A.18) in Hewitt et al (2001). 


Appendix B. The Weyl curvature tensor 


In this appendix we give an expression for the Weyl curvature parameter W in terms 
of the Hubble-normalized variables S+, E_, E 2 , Ex, A_ and A. The Weyl curvature 
parameter is defined by 


= 


77' Tpah I Tj Tjab 


6i74 


(B.I) 


where Eab and Hab are the electric and magnetic parts of the Weyl tensor, respectively 
(see WE, p 19), relative to the fluid congruence. Let and Et^p be the components 
of Eab and Elab relative to the group invariant frame introduced in appendix A. The 
Hubble-normalized counterparts of E^p and Et^p are defined according to 


^a.p — 


E, 


aP 


HaP = 


H. 




(B.2) 


p[2 ’ ' p[2 

Since Sap and Etap are symmetric and trace-free, they each have five independent 
components. We label these in analogy with the labelling of the shear variables as 
given in appendix A, i.e. 

— |('^22 + i^Ss), 


-^ 33 ), 

Sx = ^ S 23 , S 2 = Si3, S 3 = - 1 = S- 


73^23, 


(B.3) 


73 ‘ 


12 , 


with similar equations for the frame components of Etap- It then follows from (B.l)- 
(B.3) that 


= Si 


Si 


£2 


nl + n'i + nl+nl + 


Hi 


(B.4) 
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With the frame choice made in appendix A, equations (1.101) and (1.102) in WE for 
Ea/s and in conjunction with (A.1)-(A..5), lead to 


and 


£+ = E+(l + E+) - S?. - + 1^2 + 2{Nf + 3A^), 

£_ = S_(l - 2E+) + ^ + 2V3{Nf - A^), 

£x =Sx(l- 2 E+) + 8 Af_A, 

£2 = E2(1 + E+ + V3E_), 


£3 — —•\/3E2Ex , 


H+ = -3E_Af_ - aVSExA, 

H- = -3E+Af_ - 2y3E_7V_ + Ex A, 
Tfx = -(E_ +3V3E+)A-2V3ExlV_, 

n2 = Vs E,2N-, Hs = -4E2 a. 


(B.5) 


(B.6) 
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